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We discuss correlated lattice models with a time-dependent potential across a barrier and show 
how to implement a Josephson-j unction- like behavior. The pairing occurs by a correlation effect 
enhanced by the symmetry of the system. In order to produce the effect we need a mild distortion 
which causes avoided crossings in the many-body spectrum. The Josephson-like response involves 
a quasi-adiabatic evolution in the time-dependent field. Besides, we observe an inverse- Josephson 
(Shapiro) current by applying an AC bias; a supercurrent in the absence of electromotive force can 
also be excited. The qualitative arguments are supported by explicit exact solutions in prototype 
5-atom clusters with on-site repulsion. These basic units are then combined in ring-shaped systems, 
where one of the units sits at a higher potential and works as a barrier. In this case the solution is 
found by mapping the low-energy Hamiltonian into an effective anisotropic Heisenberg chain. Once 
again, we present evidence for a superconducting flux quantization, i.e. a Josephson-j unction- like 
behavior suggesting the build-up of an effective order parameter already in few-electron systems. 
Some general implications for the quantum theory of transport are also briefly discussed, stressing 
the nontrivial occurrence of asymptotic current oscillations for long times in the presence of bound 
states. 

I. INTRODUCTION 

The Josephson effect is one of the clearest fingerprints from which one can diagnose superconductivity. A DC 
electromotive force (e.m.f.) V across a thin insulating layer in a superconductor-insulator-superconductor junction 
producesJ an AC current response. The Josephson current is given by the well known expression- 

/ = /jsin(^^t + 7o^ , (1) 

where Ij and 70 are constants. In textbooks this result is explained in terms of the order parameter at the junction; 
within the insulator one has 

^{z) = ^Piefte-^' + A^ghte"'''-''^ , (2) 

where s is the barrier width, a is a characteristic inverse length and z is the normal coordinate; V'Ze/t and tpright are 
complex constants. The quantum mechanical current / is computed by the Ginzburg-Landau formula, that is, by 
using -0 as a single-particle wave function. Thus, I is proportional to sin0, where 6 is the phase difference between 
tpieft and tpright- Letting 9{t) ^ + we obtain jj). 

The electrons cannot come into play in this picture, since the order parameter is a macroscopic concept. A very 
enlightening model argument which is microscopic in spirit is the one due to Ferrel and Pranged. Let Ht denote 
the electron tunneling Hamiltonian between the left and right superconductors. The junction with N — v pairs on 
the left and v on the right (both numbers being huge) is mapped to a tight-binding chain described by an effective 
Hamiltonian Hp that describes pair hopping. For Ht — > the energy Eq of the junction does not depend on v. So, 
the state is defined by the amplitude (p^ obtained from 

Hf(Pu^Eo(P^ + F{(P^+i+(P^^i), (3) 
where the pair hopping matrix element arising to the second-order in Ht reads 



{v + l\HT\i){i\HT\iyi 



Eq — Ei 



and the sum is over the intermediate states i of energy Ei with an electron in the barrier. The eigenfunctions of Hp 
are plane waves in v space, with eigenvalues 

e{k)^Eo + 2Fcoak. (5) 
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A current 




(6) 



is associated to the group velocity in the fictitious chain and to the real current across the unbiased junction. A 
wave packet centered at k in the fictitious chain would be accelerated by 



and this inserted into © reproduces again, seen from another angle. The description by Ferrel and Prange is still 
aimed at a macroscopic system with huge pair numbers enabling us to use the concept of a group velocity. 

The question that arises is: what happens in small systems, when the fictitious chain is very short and the group 
velocity has no meaning? The Ferrel - Prange description must be modified, however pairs will continue to hop as 
bound units, and the AC response to a DC bias is a qualitative feature which looks well suited to test any model 
for superconductivity. Being a qualitative feature, it is hard to see how it can depend on size. Moreover, the order 
parameter which is well defined only in the thermodynamic limit must be supplemented by a microscopic description 
valid in finite systems; there, such concepts like long-range order do not apply but clear signatures like lack of 
dissipation must still distinguish superconducting from normal currents. The very nature of these signatures deserves 
investigation. 

In this paper we wish to look for precursors to the Josephson behavior far from the thermodynamic limit; small 
cluster studies are motivated by conceptual as well as practical reasons. Indeed, the growing technological interest in 
nano-systems justifies the identification and search. We focus on repulsive Hubbard-like models that we want to test 
for correlation-induced superconductivity; this is relevant to the search for possible non-conventional mechanisms, 
that have been considered by several workers e.g. in the context of high-Tc materials. The choice of the geometry 
has been prompted by previous studies^ in the framework of the W — theory; this develops the role of symmetry 
in inducing pairing from repulsive interactions. Besides, we also found the conditions leading to superconducting 
flux quantization in mesoscopic systems. We are therefore interested in testing this mechanism in time-dependent 
conditions. 

Below we use a Gedanken experiment for understanding the problem. Developing this idea, one meets the difficulty 
that small systems always give AC response to a DC bias for quantum mechanical reasons; that is, one finds oscillating 
currents analogous to Rabi oscillations. However, the frequency of the oscillations is proportional to the charge, so a 
sharp effect of pairing remains, i.e. supercurrents are signaled by doubled frequencies. This is a conceptual problem 
and we shall be in position to add more to try to clarify this point in the conclusions. 

There are further reasons of interest in this kind of problems. Transient current response are now actively studied, 
and within the TDLDA the asymptotic behavior of currents yields the static current-voltage characteristioi. Here we 
shall find a class of systems where, owing to correlation induced bound states, the asymptotic behavior is oscillatory. 
There is a clear need to develop new techniques, capable of tackling such situations. 

Here we outline the plan of this paper. In Section^we summarize the theory oiW = Q pairing in Hubbard clusters. 
Section IIIII presents our one-unit model. Here we discuss the simplest case, namely a 5-site Hubbard cluster with 
square symmetry. For short, we shall refer to it as the CUO4 unit; this is known to yield pairing and flux quantization, 
and here we show that it is also a prototype system for dynamic phenomena. One must perturb the square symmetry, 
in order to mimic a barrier in such a tiny system. The current excited by a constant electro-moving force (emf) 
is oscillatory with a frequency which is exactly the same as the Josephson frequency in Eq. J^l; moreover we shall 
observe Shapiro spikes (inverse Josephson effect) which represent a DC response to an AC bias at certain frequencies. 
A constant current (with some periodic ripple) ffows in the absence of an applied potential difference, which is typical 
of superconducting systems. All these effects disappear when pairing fails, e.g. if the repulsion U is removed. In 
Section llVl we shall extend the investigation of Josephson-like currents to large systems that are rings composed of any 
number of CUO4 units, up to the thermodynamic limit in principle. For small inter-unit hopping, the many-electron 
system can be exactly mapped into a Heisenberg-Ising-like chain. Systems hosting many pairs are described. An 
effective barrier is simulated by lifting the one-body levels of one of the CUO4 clusters; the superconducting ffux 
quantization is again obtained and time-dependent solutions with a constant emf are also presented. Finally, we 
present our conclusions in Section IVI 



^hk = 2eV 
dt 



(7) 
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FIG. 1: a) The undistorted 5-atom unit with 4 flux tubes b) The distorted 5-atom unit with 1 flux tube and reduced 2-3 
hopping (the barrier). 

II. W=0 PAIRING 

In the context of High-Tc superconductivity, several approaches based on the Hubbard Hamihonian have appeared 
in the literature. The quantity 

I\{N)^E(N)^E{'N -i)-2E{N -\), (8) 

where E{N) is the ground state energy of the system with N fermions is used to evaluate the effective interaction 
between particles in the interacting ground state. A negative A means that an effective attraction is developing from 
the repulsive interactions. We refer to Ref. '3 for the general theory based on the fact that highly-symmetric clusters 
possess 2-particle singlet eigenstates which do not feel the on-site repulsion W; these are called W = pairs. The 
persistence of the 'W = 0" character in the interacting case is discussed in detail there, with the flux quantization 
properties and their group theory analysis. The CUO4 cluster with Hamiltonian 

5 

(7,i—2 ia i 

is the simplest one to give the effect; here, is the creation operator for a fermion onto the Oxygen i = 2, ..,5, d 
destroys a fermion on the central site 1; besides, €p is the O energy level, tpd is the p-d hopping matrix element, and 
U is the Hubbard repulsion. In a wide range of parameters, this model yields A < for TV = 4. 
The magnetic properties of the CUO4 cluster will be discussed below in more detail. 



4 




FIG. 2: The ground state energy E{<l)) of the Cu04 with = 4 versus flux 0, in tpd units; tox = —0.05, Ep = 3.5, U — 5.3. 
Sohd: plot of 100* {E{(j}) — 7.97) for the undistorted C4„-symmetric Cu04 cluster with equal O-O bonds and 4 equal flux tubes. 
Dashed: plot of 1000 * {E{<j)) — 7.97) for the distorted cluster. Notice that we enhanced the effect 10 times in this case for 
easier comparison. Dotted: plot of 50 * {E{<j)) — 5.126) for the same input data as the dashed line but (7 = 0. In the absence 
of repulsion there is no pairing and the 4 particle system is paramagnetic. 



III. THE ONE-UNIT MODEL 

For simplicity, CUO4 clusters are natural units to build more interesting systems. We start from the simple 5-atom 
unit and then show how we can combine many units to build large systems in Sect. IIVI In order to have closed 
circuits where flux tubes can be inserted, we consider the time-dependent Hamiltonian (see Figure 1) 

Htot = Ho + Hiit), (10) 

where 

^iit)^ ^ii,t+ipl„Pi+i,o- +h.c. (11) 

Henceforth tpd will be our energy unit; we could dispense from introducing ep, but in this section ep — 3.5 to make 
the flux quantization pattern more pronounced. There are various cases to consider. 

A. Symmetric cluster 

Letting fj.i+i — tox we have the symmetric cluster of Figure 1 a). Moreover we insert flux tubes (see Figure 1 a)) 
by setting f^.i+i = ioaje*^^*'' with 7 = where ip is the flux in each tube; the flux does not affect the tpd bonds by 
symmetry. We observe a superconducting flux quantization due to the W = pairs as shown in Figure 2, solid line. A 
typical double-minimum pattern is found with sharp maxima due to the level crossing of paired states with different 
symmetry. This patter is obtained with tox — —0.05; if the absolute value of the hopping tox is increased too much, 
the coupling to the flux becomes excessive compared to the pairing energy and the system turns paramagnetiA There 
are two minima aX ip — ^ and = 0o; in both minima, there is pairing (A < 0). This pattern disappears together 
with pairing for U = 0. 

For a time-independent flux, the t bonds carry a screening current, since j oc ((> — 0(0 is equivalent to 

applying an emf across each 0-0 bond. In the absence of interactions, the current response to a constant emf across 
e.g. the 2-3 bond is periodic due to the finite system (see below). If the flux varies slowly, the system will follow the 
ground state curve adiabatically up to the crossing point. However since each level has conserved quantum numbers. 



FIG. 3: The dynamical current response in units of (solid line) versus time flowing in the 2-3 bond with (p{t) = r]t,r] = 
This value of r] is below the gap in the spectrum, which is of order \tox \ = 0.05. 



the system cannot acliabatically follow the ground state with slowly increasing (/>. As a result, the adiabatic current is 
periodic with the same period as in a normal system (e.g. the cluster with U = 0). 

B. Static properties of the mildly distorted cluster 

A more interesting behavior is obtained by breaking the square symmetry, as shown in Figure lb). A supercon- 
ducting pattern is still visible (Figure 2, dashed line). Here, = tg^ is a constant for all bonds, except for the 
bond 2-3, where we take ^2,3 = ^e^^^*\ No flux is present in the other 0-0-Cu triangles. In Figure 3, dashed line, 
wc sec that the superconducting double minimum pattern is still there but the level crossing is avoided opening a gap 
in the 4-fermion spectrum. However, the repulsion is needed to produce pairing and flux quantization. The dotted 
line shows the drastic effect of setting U = 0: the system becomes paramagnetic and normal. 



C. Emf and frequency of currents 

In the paired situation of Figure 2, dashed line, we will consider the timc-dcpcndcnt case, (p/cpo Vt, where V 
plays the role of an electromotive force, and t in units of . An adiabatic behavior is now achieved when the time 
dependence is slow compared to the gap. We are using the Gauss system, hence 

<^o = - (12) 
e 

and the Faraday law reads 
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FIG. 4: The dynamical current response in units of (solid line) versus time flowing in the 2-3 bond with </>(t) = rjt, rj 
as above, but (7 = 0. 



Therefore, we introduce the frequency rj by setting 



V = 



eV 



(14) 



Equation becomes 



Ij 



sin(2?7i) 



(15) 



for superconductor junctions, while the 2 is missing for normal systems. 

This estimate agrees with the results of Figure 3. The solid line is the current through the 2-3 bond obtained 
by integrating the Schrodinger equation. The initial conditions are = with the system in the ground state with 



iV=4. The dotted line, reported for comparison, is the adiabatic response jad{t) = c 



Og(4) 

a0(t) 



The two curves agree fairly 



well. Since tox *C 1, the main contribution to the current response comes from the 123 plaquette, which responds 
with frequency 77 for normal systems and 2irj for superconducting behavior. Some current contribution involves other 
plaquettes as well, perhaps mixing different frequencies. The agreement with the adiabatic approximation confirms 
the present analysis and will be used in Sect. IV . The normal response would be periodic with half the frequency. 

The system becomes normal for U — Q, a. fact which is apparent from the magnetic response (dotted line of Figure 
2). The current response through the 2-3 bond also become normal, as we can see from the halved frequency of the 
current response in Figure 4. Note that the adiabatic result (dashed) does not agree very well with the response 
obtained by integrating the Schrodinger equation (solid). In fact, the gap ~ tox is missing in the spectrum. 



D. Supercurrent flowing without bias 



In a macroscopic junction one can excite a supercurrent that will last forever if the exciting bias is removed, unless 
there are dissipative elements elsewhere in the circuit. This is also borne out in the paired situation of Figure 2, 
dashed line; we will consider the time-dependent case, 4>l 4>o — > Vt, where V plays the role of an emf of finite duration. 
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FIG. 5: Current response in units of of the mildly distorted Cu04 unit. A constant emf across the 2-3 bond (barrier) 
lasts from t = to t = 20. At t = 0, The current response to a flux (p{t) — jo + "n't where r; = as above and 70 = f. A 
constant supercurrent is superimposed to the oscillations. 



In Figure 5 we plot the current flowing across the 2-3 bond (barrier) when the emf lasts from < = to t = 20. At 
t — 0, (f) — ^ . A constant supercurrent is superimposed on the oscillatory response. 



E. Shapiro-like effect 



Still in the paired situation of Figure 2, dashed line, we also consider the time-dependent case of an oscillating field. 
In superconducting junctions there are spikes in the current-voltage characteristics when the voltage has a constant 
component and an oscillating one at radio- frequency lu^ of amplitude Vr (Shapiro effect): 



Vtot = V + Vr COs{LUrt). 



Then Eq.^ reads: 



/(t) =/jsin 



2eVt 



■70 



2eVr 

TlLOr 



(16) 



(17) 



In particular one finds a zero-voltage spike, because at = the resulting I{t) has a DC component. We can mimic 
such a behavior with our model. In Figure 6 we plot the current response through the 2-3 bond to a flux 



4>{t) =70 + 2 — sin(a-'rO 



(18) 



with 7o 



2' 



200 



and u)r 



10 



We observe the typical Shapiro- like DC current with a ripple on it. 



F. Comment on the one- unit model 



It appears from the above that the 5-atom cluster with suitable parametrization can simulate all the characteristic 
signatures of Josephson junctions, as if it already contained a seed of the order parameter. This is indeed surprising, 
but we show below that combining the small units we can go smoothly to large systems still keeping the same 
characteristic behavior; simplicity stems from the fact that if there are many pairs they move coherently. 
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FIG. 6: The current response in units of to a flux (f>{t) = 7o + §^ sin(cjri) where 70 



constant component is superimposed to a ripple due to the alternative paths available to the current. 



|- and tJr = Yff. The Shapiro-like 



IV. THE MANY-UNITS MODEL 



In This Section we study the time-dependent propagation of many W — pairs along rings of CUO4 Hubbard 
clusters in presence of a barrier. 

Our system is made of L — 1 CUO4 units plus a "barrier" standing between the 1-st and the L — 1-th units. The 
barrier itself is taken to be a CUO4 cluster, but with very high on-site energies, in order to make its occupancy only 
virtual. 

The L — 1 units are connected each other by the usual inter-cluster hopping Hamiltonian Ht, with real hopping 
integral r, while the hopping between the barrier and the 1-st and the L — 1-th units is assumed to be complex and 
time-dependent: r'e'^'-*-', where the phase 7(t) will be discussed later. By means of such a choice we have in mind 
that in transport gedankenexperiments the potential drops only through the barrier. 

According to the notation of our previous work^, the model reads: 



Hq + + H] 



barrier 



(19) 



with 



L-l 



(d) (d) , Y^^(p) ^(p) 



where, p]^ is the creation operator for a hole onto the Oxygen i = 1, .., 4 of the a-th cell and so on. In this section 



ep = 0. 



Hr is an inter-cell hopping Hamiltonian which allows a fermion in the i-th Oxygen site of the a-th unit to move 
towards the i-th Oxygen site of the a ± 1-th unit (without involving the L-th unit, which is the barrier) with hopping 
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integral t: 

L-2 
a — 1 ia 

Finally we model the barrier with the following Hamiltonian: 

where Ntot is tlie number of particles in the L-th cluster and W is an extra (large) energy felt by each particle in 
the barrier, which can be provided, for instance, by a gate voltage. The barrier is linked to the rest of the ring by a 
complex hopping Hamiltonian: 

Hr,{t) = r'e^-'^*) Y.^pl_,,,PL,^a +pU,Pi.a) + h.C. = T'e'^^'^'Hj + h.C. (21) 
In the static case we set (with the same notation o^) 

27r0 , , 

A0O 

where is the total flux through the ring and A = 2 as long as the flux is associated with two bonds. In the 
time-dependent case, we will set (/i/^o — Vt, where V plays the role of an electromotive force, as in the previous 
Section. 



A. Static Case 



Let us begin with the static case. As usual we introduce + Ht-i perturbatively, in the same spirit of our previuos 
workS. When U/t :^ 5, t = t' = and the number of particles is appropriate and even, each CUO4 cluster is populated 
by 2 or 4 particles (since the occupations with 3 particles have a gap of the order of |A| ), while the barrier is totally 
empty, since the extra energy even for a single occupancy is very large. Now, if |r|, |r'| << |A|, the perturbations 
are operative only by the second order. In particular Hr provides an effective hopping r^ff t^/|A| for bound 
W = pairs, which can propagate in the l,2,...,i — 1 clusters as hard-core bosons. On the other hand Hr' acts 
differently. Indeed the occupancy of barrier is highly unfavoured in the ground state because of the large energy W. 
As a consequence the W — pair can jump from the 1-st cluster to the L — 1-th cluster an viceversa only at the 

fourth order in Hr' . In this case the effective hopping through the barrier is Je *o with 

(23) 

where 

e = E{3)+E{2)+E{l) + W (24) 
(with E{N) the ground state energy of the CUO4 with N particles), and C is a order 1 costant given by 



C = (*L-i(4)| ® (*i(0)| (*i(2)I F7|*i(2)) ® |*i(l)) <S> |*L-i(3))x 
(*L-i(3)| ® ® (*i(2)| i?j|*i(3)) |*i(0)) ® |*i(3))x 

(*L-i(3)| ® (*l(0)| ® (*i(3)| ff/|*i(3)) ^ |*l(1)) ® |*L-i(2))x ^^^^ 
(*L_i(2)| ® (*i(l)| ® (*i(3)| -ffj|*i(2)) ® |*l(0)) ® |*i(4)), 

where \'i>a{N)) is the ground state of the a-th cluster with N particles. 

The resulting low-energy effective hamiltonian is then equivalent to a 1-dimensional spin 1/2 chain with C = L — 1 
sites, given by (see RefiiSi for the case of r' = 0): 

£-1 £-1 

H^fficf,) = -2T,ff ^ 5^5^+1 + T,ffY.{S+S-+, + h.C.) - T^nSlSl + (Je'^5+5r + h.c). (26) 
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FIG. 7: E{(j>) — -E(O) (in units of tpd) as a function of <j}/<j}o for 2 x Sfot ~ Ij 2, 3, 4, 5, 6. Here we have taken and J/r = 0.1. 



We observe that the barrier has been "integrated out" and the resulting model has the 1-st and the C-th sites linked 
by J. 

It is worth to note that this model, because of the the presence of two different parametres Tg// and J, is not 
solvable by the Bethe ansatz, therefore we must resort numerical methods. We solved Heff{4>) with C up to 12 and 
Sfgf. up to 3. In this case, in the original fcrmionic model there were L = 13 CUO4 clusters (12 -I- barrier) and 6 
W — pairs, that is 12 x 2 + 6 x 2 = 36 fcrmionic particles. In other words the number oi W = pairs is 2 x S'fgj, 
which is a conserved quantity. 

In figure we show the ground state energy E{(j)) as a function of the magnetic flux for C = 12 and 2 x S'f^j — 
1,2,3,4,5,6. It is remarkable that whatever is the number of = pairs in the original systems (1,2,3,4,5,6), the 
peridicity of E{(j)) does not change and shows the same superconducting flux quantization. This is a very interesting 
feature, suggesting that all the W = pairs may behave " coherently" , providing the same macroscopic response, no 
matter of their number. It is also clear that the oscillations of E{(j)) get more and more pronounced as the number of 
pairs in the ground state become of the same order of number of the lattice sites. This is an indication that such a 
behavior should survive in the thermodynamic limit. 

Finally we observe that the trend of E((j)) by incrcsing the number of particles is deeply different with respect to 
the one observed in the repulsive Hubbard ring. In the latter case^"'^^, there are many minima between (j) — and 
(/> = 00, depending on the ratio Nt/LU, where N is the number of particles and L is the number of lattice sites. 
On the other hand, in the case under consideration it seems that the superconducting nature of the charge carriers 
"freezes" the period oscillations to be (j)o/2. 

As a last remark, we notice that E{(j)) is a smooth function of (p. This means that no level crossing occurs between 
(p = and (j) = 4>o- This is a direct consequence of the breaking of the translational symmetry induced by the presence 
of the barrier. As discussed in the previous Section, this condition will be important in the time-dependent case, 
where, at least in the adiabatic approximation, the system at t = (supposed to be in the ground state at = 0) will 
be able to evolve continuosly to the static ground state at = 0o/2. This condition is necessary in order to have a 
time-dependent superconducting response. 



B. Time-Dependent Case 

Now we consider a time-dependent perturbation, with the effective hopping through the barrier given by Je*'*'^^*. 
In this case we rely on the adiabatic approximation, such that, nevertheless of the time-dependence, the low energy 
Hamiltonian is the same as Heff, with (j)/(f>o — * Vt- The correctness of this approximation has been checked directly 
in the previous Section, in the time-dependent analysis of the single CUO4 cluster. There, the exact time-dependent 
solution of the whole original electronic problem was achieved and compared to the approximated one. 

The exact time-dependent solution of the Schrodinger equation associated to Tfe/ / (t) is much harder than the static 
solution. Therefore we only consider the case C — A and 2 x S^^^ — 2. Anyway, in the original fcrmionic model, it 
means L = 5 CUO4 clusters (4-|-barrier) and 2 = pairs, that is 12 particles. As initial condition at t = we take 
the gruond state of Hf,ff{4>) at (/> = 0. 
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FIG. 8: I{t) as a function of Vt for 2 x Stot = 2 for two intervals (0,1) and (49,40): Here we have taken and J/r = 0.1. 

Our aim is to evaluate the time-dependent current I{t) through the barrier, i.e. through the bond governed by 
Je *o J according to 

lit) = J(*(i)|e'4"^*5+5f - e-'^^'^'S^S+l^it)). (27) 

The exact numerical solution shows that I{t) is a peridic function of t with period 1/{2V). In figure we plot I{t) 
as a function of Vt in two different unitary intervals: (0,1) and (48,49). 

We stress that the time-dependent response of the many-units model have the same Josephson-like behaviour of 
the single-unit model. 

V. CONCLUSIONS 

We have shown how to model a Josephson-like junction using microscopic ingredients like electrons hopping in an 
atomic lattice, a tunneling barrier, and on-site correlation. The barrier consists of a bond characterized by a weak 
time-dependent hopping, where the emf is applied. The current response is easier to calculate and analyze when 
there is emf on a single bond. We stress that the Hubbard interaction U leads to a current response in quantitative 
agreement with Eq. This happens in circumstances that we understand, namely: i) the formation of W = 

pairs, ii) a well developed double minimum in the energy-flux plot, iii) a soft distortion of the symmetry opening 
a small gap (due to avoided crossing) in the many-body spectrum, iv) a small enough emf V, in order to allow for 
adiabatic response. When the last condition is realized, calculations based on the adiabatic approximation are in 
good agreement with the full solutions; this allows us to extend the analysis to fairly large systems. In modeling rings 
of Hubbard clusters, the barrier was formed by a Cu04 unit with high one-electron energy; in order to deal with the 
low-energy physics, we derived an effective Hamiltonian where the barrier is integrated out and represented by an 
effective, weakened, complex bond, similar to the barrier in the above distorted CUO4 case. 

In the single-unit model, the superconducting AC current has the same frequency 2^ as one would expect to 
observe in a macroscopic junction; in the latter case, however, the normal response is a DC current. In our cluster 
calculations, on the other hand, the current response of the normal systems is periodic with half the Josephson 
frequency; this happens with one-particle systems and with 4-body ones with U = 0, when no pairing occurs. This is 
due to a size effect that we wish to comment upon here. 

As shown in Ref. l)llh . in the non-interacting case (which is clearly normal) the transient current across a tunnel 
barrier between infinite leads tends to a constant for long times; the asymptotic t ~> 00 behavior yields the current- 
voltage characteristics calculated by several people by time-independent approachesiSMMi^. This has recently been 
proven as a theorem by Stefanucci and Almbladh^ within the Time-Dependent Density Functional Theory in the 
Local Density Approximation. In order to understand the difference, we recall that the infinite junctions considered 
in (— ) never reach equilibrium, and the steady current that holds in the t — > 00 limit is just a never-ending transient. 
On the contrary, the oscillating behavior of the small cluster follows almost adiabatically the instantaneous ground 
state. Another difference is that the normal system is unable to trap the magnetic ffux: when the field penetrates 
the sample our modeling should be less appropriate than in superconductors. So we can expect that the normal 
current oscillation at half frequency, that we observe in the small cluster, should acquire a longer and longer transient 
behaviour with increasing sample size, while the field penetration should superimpose different oscillation frequencies. 
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Instead, the superconducting pairs should keep their frequency in Hne with equation Q right to the thermodynamic 
hmit. 

According to our interpretation, the Josephson effect is essentially another facet of the superconducting flux quanti- 
zation. The necessary distortion is always present in macroscopic systems, where in addition a large energy separates 
the states of different pair symmetry. Unbiased supercurrents and Shapiro spikes are also obtained by the same models 
that mimic the AC Josephson behaviour. 

In the many-units model, the CUO4 clusters are combined in large systems linked by weak hopping integrals; the 
solution is found by mapping the problem into a spin 1/2 chain. In the time- independent case, we presented evidence 
that the flux quantization pattern is unmodified when several bound pairs propagate together, except that the energy 
barriers separating the minima grow large. In the time-dependent case we solved Schrodinger equation numerically, 
and again we got a Josephson-junction-like behavior. 

Finally, it is remarkable that even with one pair, that can be on the left or on the right, there is enough uncertainty 
on the pair number to create some ghost order parameter. 
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